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Q^. ABSTRACT 

6 

H 

'o^ I We present Hobby-Eberly Telescope (HET) observations for galaxies at red- 

' shift z < 0.3 from the Sloan Digital Sky Survey (SDSS) showing large velocity 

dispersions while appearing to be single galaxies in HST images. The high signal- 

>• ' to-noise HET spectra provide more definitive velocity dispersions. The maximum 

velocity dispersion we find is a* = 444 km s~^. Emission-line widths in QSOs 



O ' indicate that black holes can exist with masses M, exceeding 5 billion Mq, im- 

• ■ plying cr* > 500 km s~^ by the local M, — cr* relationship. This suggests either 

O ■ that QSO black hole masses are overestimated or that the black hole - bulge 



relationship changes at high black hole mass. The latter option is consistent 
with evidence that the increase in a* with luminosity levels off for the brightest 
elliptical galaxies. 

Subject headings: galaxies: general — black hole physics 



1. Introduction 

The use of broad emission lines in QSOs to measure black hole mass, M,, allows study 
of black hole demographics over a large volume of space. The largest black hole masses range 
up to 10^° Mq, exceeding the largest M, {10^-^Mq) measured for nearby galaxies. Shields et 
al. (2006a) derive a space density of 200 Gpc~^ for black holes with M, > 10^'^ Mq based on 
broad line widths in the most luminous QSOs. Wyithe & Loeb (2003) find a similar result 
based on the assumption that all QSOs are radiating at their Eddington limit. However, from 
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the M, — (J* relationship (Gebhardt et al. 2000a; Ferrarese & Merritt 2000; Treinaine et al. 
2002), a black hole with mass 10^''^ Mq corresponds to a velocity dispersion of 500 km s~^, 
larger than any a^ observed in local galaxies. Indeed, in a volume of space ~ 0.5 Gpc'^ 
(approximately corresponding to the SDSS DR2; Strauss et al. 2002) none appear to have 
(j*>430 km s~^ (Bernardi et al. 2006, 2008, hereafter BOG and BOS respectively). These 
seemingly contradictory results raise the question of where the largest black holes reside, 
and what value of a* characterize their host galaxies. 

Accurate knowledge of the largest values of a* that occur in galaxies is therefore im- 
portant for the understanding of the largest black holes in the universe, as well as having 
obvious bearing on the evolution of the largest galaxies. However, the signal-to-noise (S/N) 
of survey spectra — especially in the case of higher redshift galaxies — is often too low to 
make an accurate determination of cr* or to confidently rule out doubles and superposition 
as causes for large apparent velocity dispersions. With this goal in mind, we have selected 
eight galaxies which have amongst the highest velocity dispersions in the sample of BOS and 
appear, in i/ST images and SDSS spectra, to be single galaxies. We have obtained high S/N 
spectra with the Hobby-Eberly Telescope (HET) to measure a* with the greater accuracy 
and to double-check the likelihood of superposition. 

Absolute magnitudes used in this study are calculated using the cosmological parameters 
i/o = 70 km s"^ Mpc-\ ^m = 0.3, and ^a = 0.7. 



2. Sample Selection 

B06 have studied high-a* galaxies in a sample from the SDSS approximately correspond- 
ing to Data Release 2 with z < 0.3, corresponding to a co-moving volume of approximately 
0.5 Gpc~'^. The selection criteria are described in detail in B06. Briefly, they chose early-type 
galaxies with Petrosian apparent magnitudes in the range 14.5 < rpet < 17.75 and measured 
velocity dispersions in the SDSS spectroscopic pipeline. Out of a total of ~ 40, 000 galax- 
ies, they found 100 candidate galaxies with a* > 350 km s~^. Of these, roughly half are 
superpositions as shown by the SDSS images or spectrum line profiles. For most of the 
remaining objects, BOS obtained HST images that in turn reveal roughly half of them to be 
superpositions. The surviving 23 objects have values of o"^, ranging up to ~ 430 km s~^ as 
measured from the SDSS spectra. We selected eight of these objects based on their large 
(J* and availability in the sky during the proposed observing period. Figure [1] shows HST 
images of these galaxies from BOS. 



-3- 



3. Observations and Data Reduction 

Observations were obtained using tlie Marcario Low Resolution Spectrograpli (LRS; 
Hill et al. 1998) on the 9.2-ni Hobby-Eberly Telescope during the period spanning April - 
November, 2006. We used a slit width of 2", a 600 line mni"^ grating, and a GG385 blocking 
filter; the 3072 x 1024 Ford Aerospace CCD has a plate scale of 0.235 arcsec/pixel and was 
binned 2x2. Exposure times were 40 minutes each. Wavelength coverage for this setup 
is 4300 - 7300 A, with a resolving power R = A/AA = 650. The instrumental width is 
o'inst = FWHM/2.35 = 200 km s~^. On most nights the sky conditions were spectroscopic 
with 1.5 — 2.0 arcsecond seeing. Table [2] lists the spectral extraction apertures and galaxy 
morphologies. 

The data reductions were carried out with FORTRAN algorithms developed by one of the 
co-authors (K.G.) using standard reduction techniques. We first performed CCD corrections 
for overscan. Flat fields were created from five images of an internal continuum source and 
normalized, and each image was then divided by the appropriate fiat field. Rectification 
was performed by determining the x and y locations of the spectra on the images and 
approximating a trace. Images were sky subtracted using apertures for the night sky that 
were defined manually for each object. Images for each galaxy were combined and then a 
1-dimensional spectrum was extracted by defining an aperture for the object and summing 
pixel values for each wavelength bin. The wavelength scale was calibrated by means of a 
fourth-order polynomial wavelength solution for an emission spectrum of a Cd lamp. 



3.1. Measuring a^ 

Details of the spectrum modeling are described in Gebhardt et al. (2000b). Briefly, 
we used an automated procedure to fit the spectrum separately in the regions of the G- 
band line and the Ca II H+K lines using an adjustable combination of stellar templates 
to simulate the spectrum of an elliptical galaxy. Table [1] lists the template stars and their 
spectral types. There is some difficulty in simulating spectra for large-a* galaxies, because 
the stellar population is not well known. For this reason, we constrained our fitting routine 
to the regions immediately surrounding the G-band and Ca II H+K lines. A measurement 
for cr* was obtained through a simultaneous fitting of the continuum and the velocity profiles 
of the G-band line or the Ca II H+K lines until an optimal match to the galaxy spectrum 
was achieved (see Figure [2]). We modeled the velocity proffies using two methods: (1) 
least square fits of a pure Gauss function; and (2) least square fits of a Gauss-Hermite 
function, which is Gaussian multiplied by Hermite polynomials. The coefficients of the 
Hermite polynomials, h^ and h^, characterize the asymmetrical and symmetrical deviations 
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from a pure Gaussian profile, respectively (see van der Marel and Franx 1993). The first 
two moments of the Gauss- Hermite represent the mean velocity and velocity dispersion, a*. 
We quote the a* measurements from fitting method (1), and used the parameters of the fit 
from fitting method (2) as a check for binarity or other irregularities. Goodness of fit was 
confirmed by visual inspection of the model fit and velocity profiles for each galaxy. 

Error bars were determined via Monte Carlo trials in which the initial data were ran- 
domly varied based on the rms noise per pixel and then re-fitted by our automated procedure 
to get a 0"=,,. The standard deviation of a,, from 100 such trials determined the la level of 
error. The errors for the G-band fits and Ca II H+K fits are comparable, but we quote only 
the cr* of the G-band measurement, since it is less susceptible to template mismatch due to 
a-element enhancement (see Barth et al. 2002). 



4. Results 

Table[3]lists a* measured for each galaxy from both the HET and SDSS spectrum. Figure 
|3]shows that in most cases there is agreement within ±10%. The largest a* we find is 444 ±16 
km s~^ for SDSS J082646. 72+49521 1.5, which differs from the SDSS spectrum measurement 
by ~ 9%. The smallest a, in our sample is 302 ±18 km s~^ for SDSS J082216.57±481519.1, 
which differs from the SDSS spectrum measurement by ~ 25%. Figure [3] shows that the 
largest discrepancies are for lower-cr* objects. We apphed our fitting algorithm to the SDSS 
spectra, and our measured dispersions agree with those reported in B06. There is no obvious 
reason for the discrepancies. One possibility may be that the galaxies were selected for large 
cr* from SDSS. Given the steeply dropping number of galaxies with increasing a*, and the 
large number of galaxies in the full SDSS dataset, noise may act to elevate a few objects 
from 0"* ~ 300 km s~^ to a measured value of a* >400 km s~^. 



4.1. Tests for Binarity 

We attempted to assess the likelihood of binarity by inspecting plots of the velocity 
profiles for asymmetries and multiple components — fiat-topped profiles in particular indicate 
binaries. The strength of the Hermite coefficients h^ and h^ served as numerical indicators 
of binaries: non-zero h^ indicates asymmetry, with /13 < corresponding to a blue wing on 
the profile and h-s > corresponding to a red wing; h^ < indicates a boxy or fiat-topped 
profile. The {h^l and {h^l for each of the galaxies in our sample were modest (< 0.1), and none 
of the velocity profiles showed indication of asymmetry or multiple components. However, 
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this method of assessing binarity appears to detect only obvious binaries with very large 
measured o"^. We fit the SDSS spectra of the galaxies in Table 2 of B06, which are flagged 
as superpositions, to see whether h^, and /i4 would help to distinguish them as binaries. We 
found that only in cases where measured a* >500 km s~^ and visual inspection of the profile 
fits showed obvious signs of binarity did the \h^\ and \hi\ exceed 0.1. We rely on the HST 
images and lack of gross effects on the line profiles to rule out binaries, but cannot rule out 
precisely superimposed images and modest velocity separations that could enhance measured 
cr*. 



4.2. The a^-L Relationship 

In Figure HI we plot a* against V-band absolute magnitude and compare our galaxies 
to the cr* — L relationships and data from Figure 3 of Lauer et al. (2007a), which includes 
galaxies with core, power-law, and intermediate type profiles. Lauer et al. point out that 
the absolute magnitudes given for these galaxies by the SDSS catalog are underestimated 
due to improper sky subtraction in the SDSS pipeline. Hyde et al. (2008) corrected for 
this by carrying out their own photometric reductions of these objects, so we follow BOS 
in adopting their r-band absolute magnitudes, and convert to Johnson V-band magnitudes 
using the prescription from Jester et al. (2005) iV = g — 0.59[5' — r] — 0.01). The four lowest 
luminosity galaxies agree with the a^ — L relationship for power- law galaxies (L ~ cr^) , while 
the higher luminosity galaxies agree with the relationship for core galaxies {L ~ al) — a direct 
analysis of the HST light profiles confirms that the low L objects are power-law galaxies, 
and the others are core galaxies (Hyde et al. 2008). Figure [5] shows the residual from the 
Faber- Jackson (L ~ aj) relationship versus the residual from the effective radius-luminosity 
relationship (Bernardi et al. 2007a). More concentrated bulges are correlated with greater 
deviation from the Faber- Jackson relationship, in the sense that our galaxies have measured 
velocity dispersions larger than what is predicted from their luminosities. See Section ^ for 
the significance of this correlation. 



4.3. The M, — a^ Relationship 

Table [3] also lists black hole masses predicted from the M, — a* relationship for nearby 
galaxies quantified by Tremaine et al. (2002) and from the log-quadratic formulation of 
Wyithe et al. (2006). For our highest-a* object, we predict M, k. 10^'^ Mq (Tremaine) 
and K. 10^-^ Mq (Wyithe), rivaling M, for M87 (Harms et al. 1994) and recently that of 
Abell 1836-BCG, the most massive black holes measured locally (Dalla Bonta et al. 2007). 
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Note that the a* of 375 km s~^ for M87 (Tremaine et al. 2002) places it above the M, — a^ 
relationship by a factor of about two in M,. If the M, — a* relationship changes for high 
M, (e.g. Netzer et al. 2003), black hole masses could be larger for our galaxies than what 
is predicted here. Alternatively, M, may not correlate well with any galaxy properties at 
high M,. Given the scatter in the black hole - galaxy relationships, and lacking direct 
measurement of black hole masses > 10^'^ M©, the nature of the M, — a* relationship 
remains uncertain for large M,. 



5. Discussion 

The largest a^, we measure is ~ 444 km s~^ in a volume of space corresponding to 
0.5 Gpc^. We find no observational support for galaxies with a* > 500 km s~^, contrary 
to the QSO-based prediction of 100 such galaxies (see ^ID- Evidently, either (1) black hole 
masses derived from QSOs are overestimated or (2) the linear relationship between log M, 
and log cr* breaks down for high M, (Netzer 2003). Regarding the first possibility, QSO 
black hole masses are typically estimated using the "photoionization method" (see Shields 
et al. 2003 and references therein), in which the virial mass is given by M, = v^RjG. The 
velocity is given by the width of broad line region (BLR) emission lines such as H/5, and 
the QSO luminosity provides an estimate of the radius of the BLR. Assuming the BLR gas 
is indeed virialized, the main uncertainties are the BLR radius-luminosity relationship and 
the geometry of the BLR (flattened or spherical). The -Rblr — ^ relationship was calibrated 
using the reverberation mapping technique on a sample of lower-luminosity AGN (see review 
by Kaspi 2007) and scales as -Rblr ~ ^'^ with 7 = 0.5 — 0.7. Bentz et al. (2006) find 
7 = 0.52 after correcting L for host galaxy star light, consistent with photoionization physics 
(Shields et al. 2003). This relationship holds for a relatively wide range of luminosities up to 
-^boi = 10^^ erg s~^ , and is extrapolated to higher luminosities to estimate black hole masses 
for bright QSOs. As discussed by Salviander et al. (2007), the calibration used by Shields 
et al. (2006a) is consistent with observed M, and a^ in Seyfert galaxies (Onken et al. 2004). 
The calibration recommended by Bentz et al. (2006) gives M, 0.1 dex larger than that of 
Shields et al. (2006a). Shields et al. point out that a majority of the highest luminosity 
QSOs in the sample of Shields et al. (2003) have M, > lO^'^M©, which argues against scatter 
in the derivation of M, from AGN luminosity and broad line width as a major source of the 
large M, for this sample. Thus, an overestimation of M, in QSOs appears unlikely as a 
cause of the conflict between the abundance of large black holes in QSOs and the number of 
galaxies with sufficiently high a*. 

We must therefore consider the possibility of a breakdown in the linear relationship 
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between log M, and log a* at high M,. This is consistent with the properties of the brightest 
cluster galaxies (BCGs). As shown in Figure H] and discussed by Lauer et al. (2007a), the 
Faber- Jackson relationship levels off for the most luminous galaxies (also see Oegerle & 
Hoessel 1991; Bernardi et al. 2007a). This suggests that M, may track galaxy luminosity 
rather than a* for the largest galaxies. Lauer et al. (2007a) have augmented the galaxy 
luminosity function of Blanton et al. (2003) to include brightest cluster galaxies (BCGs) 
from the Postman & Lauer (1995) sample. This luminosity function (Lauer et al 2007b, 
Figure 1) gives 10^'^ galaxies per Gpc^ with My < —23.9. If current black hole samples 
are unbiased (see discussion in Bernardi et al. 2007b), then this luminosity corresponds to 
M, = 10^'^ Mq (see Equation 6 of Lauer et al. 2007a). Within the uncertainties, this is 
an adequate number of BCGs to host the largest black holes observed in QSOs. Note that 
Tundo et al. (2007) show that scatter in the M, — Lbuige relationship can bias estimates 
like this upwards by factors of several. A consequence is that the largest black holes will 
typically reside in proportionally modest-sized galaxies (see discussion in ^. This may 
introduce scatter into the black hole - host galaxy relationships for the largest black holes. 

Does our highest velocity dispersion represent a physical limit for 0"=,,? The a* — L 
relationship flattens out for a^ at high L, which may be explained by influences on cr* in 
galaxy merger models; the largest galaxies are believed to have formed through a succession 
of dissipationless, or "dry," mergers. Loeb & Peebles (2003) discuss possible causes of a 
limiting value of a*. Major merger simulations by Boylan-Kolchin et al. (2006) show that 
the type of trajectory during galaxy mergers determines the concentration, and therefore a*, 
of the remnant spheroid. For orbits with high angular momentum, more energy is transferred 
from the stellar component to the dark matter halo, resulting in a more concentrated remnant 
and a larger a*. Figure [5] is consistent with this picture, showing more concentrated bulges 
(in the sense that the measured effective radius. Re, is smaller than that predicted from 
the Re — L relationship) corresponding to over-large a* (for a given bulge luminosity and 
the L ~ cr^ relationship). Since we use luminosity to predict both a^ and Re, the axes in 
Figure [5] are not independent. However, these correlated errors cannot explain the trend in 
the figure. Given the definition of the residuals, A log cr* = log cr*(L) — log cr* and A log 
Re = log Re{L) — log Re, and that Re oc L^'^ and cr* oc L°'^^, errors in luminosity will 
affect A log Re and A log cr* in the sense of creating a positive slope in the figure rather 
than the negative slope that is actually seen. The magnitude and sense of the change in the 
relationship due to correlated errors is indicated in the figure. The negative slope in Figure 
[5] is qualitatively consistent with expectations from the virial theorem — for a given spheroid 
mass (and therefore luminosity) a decrease in radius corresponds to an increase in measured 



6. Comment on Evolution of the M, — cr* Relationship 

Recent studies by Shields et al. (2006b), Peng et al. (2006), Salviander et al. (2007), 
and others suggest that the M, — a* relationship evolves with redshift in the sense that, for 
a given M,, the velocity dispersion is too small at higher redshifts. Whether this apparent 
evolution is real is controversial. For example, Salviander et al. (2007) studied M, and a* 
in SDSS QSOs at redshifts up to 2 = 1.2, estimating a^, from the width of [O ill] A5007 and 
[O ll] A3727. They found that higher redshift QSOs had M, too large for a* by ~ 0.5 dex. 
The interpretation of this trend requires caution. Brighter QSOs generally have larger black 
holes, so that larger M, corresponds in practice to higher redshift and larger look-back time. 
(1) Salviander et al. (2007) and Lauer et al. (2007b) describe how a Malmquist-like bias 
leads to over-representation of large black holes in modest galaxies. This bias is greater at 
the steep, high luminosity end of the QSO luminosity function; and it may give rise to the 
appearance of evolution in the M, — a* relationship. (2) An observational bias favoring the 
detection of narrow lines affects the statistics of the [O ill] and [O ll] line widths. Salviander 
et al. find that these two biases account for about 0.3 dex of their apparent evolution, leaving 
~ 0.15 dex of residual evolution. However, as discussed above, a* levels off in the cr* — L 
(Faber- Jackson) plot (Figure Hj) at high luminosity, where brightest cluster galaxies (BCGs) 
do not follow the relationship (Oegerle & Hoessel 1991; Lauer et al. 2007a; Bernardi et al. 
2007a). Use of the local M, — a* relationship (Tremaine et al. 2002) with M, oc cr^ at high 
M, could overestimate the expected a^, and give the appearance of evolution. A correction of 
only 0.1 dex on this basis applied to the residual evolution of Salviander et al. (2007) would 
leave no significant remaining evolution. Alternatively, Bernardi et al. (2007b) suggest that 
there is little bias in the M, — o"* relationship. Instead, they suggest that bias occurs in the 
M, — Lbuigc relationship, as a result of a steeper slope in the a^, — L relationship for black 
hole samples compared to the SDSS sample. 



7. Conclusion 

We have used the HET LRS to obtain high signal-to-noise spectra for eight galaxies 
selected from the sample of B06-B08 representing galaxies with the highest measured values 
for a^, that appear to be single galaxies in HST images and SDSS spectra. The maximum 
cr* we find is 444 km/s, similar to B06. We find no support for cr* > 500 km s~^. If black 
hole masses are reliable, this is indicative of a M, — cr* relationship that deviates at high 
M, from the locally-observed log-linear relationship. Though our data do not address this 
issue directly, it is consistent with the change in the relationship between a* and luminosity 
for the brightest elliptical galaxies. 
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Table 1. Template Stars 



Star 


Spectral Type 


HD10761 


G8III 


HD111812 


GO III 


HD161797 


G5IV 


HD199960 


Gl V 


HD219449 


KOIII 


HD220954 


Kl III 


HD39283 


A2 V 


HD81146 


K2III 


HD85235 


A3 IV 


HD92588 


Kl IV 


G166-45 


A5 V 



Note. — Stars contribut- 
ing to the stellar template 
used by the fitting algorithm 
to simulate the continuum of 
an elliptical galaxy. 



Table 2. Spectral extraction apertures and galaxy profile types 



Galaxy aperture aperture profile 

(SDSS J) (arcseconds) (kpc) 



082216.5+481519.1 


2.35 


5.34 


power-law 


082646.7+495211.5 


2.82 


7.74 


power-law 


124609.4+515021.6 


1.88 


7.76 


core 


133724.7+033656.5 


2.35 


5.55 


power-law 


145506.8+615809.7 


2.82 


11.8 


core 


171328.4+274336.6 


2.12 


9.38 


core 


211019.2+095047.1 


1.65 


6.08 


core 


221414.3+131703.7 


2.82 


7.51 


power-law 



Note. — The HET spectra were extracted from the apertures 
shown in columns 2 and 3. Galaxy profile types are from Bernard! 
et al. (2008). 
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Table 3. Physical parameters for galaxies 



Galaxy 


2 


Mv 


0-, 


0-, 


log M. 


log M. 


log M. 


log Re 


(SDSS J) 






(HET, G-band) 


(SDSS) 


K) 


(it*, quad.) 


(L) 








(mag) 


(km s-i) 


(km s-i) 


(Af©) 


[Mq) 


(Mq) 


(kpc) 


082216.5+481519.1 


0.127 


-21.27 ± 0.042 


302 ± 18 


402 ± 28 


8.85 


8.84 


8.28 


0.337 


082646.7+495211.5 


0.159 


-22.22 ± 0.009 


444 ± 16 


408 ± 26 


9.52 


9.69 


8.79 


0.547 


124609.4+515021.6 


0.269 


-23.95 ± 0.052 


404 ± 22 


402 ± 35 


9.36 


9.47 


9.70 


1.249 


133724.7+033656.5 


0.133 


-22.74 ± 0.004 


419 ± 29 


422 ± 31 


9.43 


9.57 


9.06 


0.645 


145506.8+615809.7 


0.274 


-24.37 ± 0.015 


408 ± 14 


394 ± 36 


9.31 


9.42 


9.92 


1.460 


171328.4+274336.6 


0.297 


-24.00 ± 0.010 


378 ± 13 


413 ± 27 


9.40 


9.52 


9.73 


1.108 


211019.2+095047.1 


0.230 


-24.37 ± 0.004 


309 ± 34 


386 ± 32 


9.28 


9.37 


9.92 


1.093 


221414.3+131703.7 


0.153 


-22.02 ± 0.008 


342 ± 18 


384 ± 28 


9.27 


9.36 


8.68 


0.391 



Note. — My calculated using r- and g-band magnitudes from Hyde et al. (2008) and the transformation from Jester et al. 
(2005). Errors for My were calculated by adding in quadrature weighted rms errors for Mr and Mg. ct* sdss source: Bernardi et 
al. (2006). M, (cr») and M, ((T*-quad.) inferred from the M, ~ o-» relationship quantified by Tremaine et al. (2002) and Wyithe 
(2006), respectively; M,{L) inferred from the M, — L relationship quantified by Lauer et al. (2007a). Re is from Hyde et al. 
(2008). 
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Fig. 1. — i^iST images from BOS for the galaxies in Table [31 Red circles show the size of the 
3" SDSS fiber. 
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Fig. 2. — Fits to the G-band region of the spectrum for the galaxies in Table [3l The left 
panel shows fits for power-law galaxies; the right panel shows fits for core galaxies. 
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Fig. 3. — cr*(SDSS) versus cr*(HET). The error bars show the la level of error. The solid 
line shows the 1:1 relationship. 



-17- 




2.5 - 






b 



.5 



■20 



'22 



-24 



-26 



M, 



Fig. 4. — Measured a* versus l^-band absolute magnitude. Large circles show data for our 
galaxies. The leftmost four circles are power-law galaxies; the remaining circles indicate 
core galaxies. Error bars show the la errors in a* and My, however error bars for My 
are smaller than the data points. The remaining data are from Figure 3 of Lauer et al. 
(2007a): squares show data for power-law galaxies; pentagons show data for intermediate 
type galaxies; triangles show data for core galaxies. The dashed line corresponds to L oc a^ 
and the solid line to L oc cr^. 
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Fig. 5. — Residuals from the bulge velocity dispersion-luminosity (Faber- Jackson) relation- 
ship versus residuals from the effective radius-luminosity relationship. Residual a,, = log 
cr*(L) — log cr,,(HET) and residual R^ = log Re{L) — log Re- Log a^{L) was calculated using 
the a^, — L relationship (Equation 2 of Bernardi 2007) and log Re{L) was calculated using 
the Re — L relationship (Equation 2 of Bernardi et al. 2007a). The arrows in the upper right 
of the figure indicate the magnitude and sense of the change in the relationship that would 
be due to changes in the predicted Re and cr* should My change by a factor equal to the 
average error in My, 5 My = 0.018. 



